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Abstract 

The CLEO Collaboration has observed the first hadronic transition among bottomonium (bb) 
states other than the dipion transitions among vector states, T(nS)^7r7rT(mS). In our study of 
T(3S) decays, we find a significant signal for T(3S)^7a;T(lS) that is consistent with radiative 
decays T(3S)-^7Xfei,2(2-P), followed by X6i,2(2-P)^Ci^T(lS). The branching ratios we obtain are 
Bixbi{2P)^i.'T{lS)) = 1.63t°;110;i6% and e(x62(2P)^u;T(lS)) = l.lOt^oM~oil%^ in which the 
first error is statistical and the second is systematic. 
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FIG. 1: The spectrum of bottomonium states. Dilepton decays of the vector states are denoted 
by the solid lines, while dipion transitions are denoted by dashed lines. The radiative decay 
T(3S)— >X6(2i^) is indicated by the dot-dashed line, and the decay Xfe(2P)^wT(lS) by the dotted 
line. 



The only hadronic decays of bottomonia that have been experimentally observed to date 
are the tttt (tttt = tt+tt^ and 7r°7r°) transitions among vector bottomonium states. 0| Fig- 
ure n illustrates the bottomonium spectrum and transitions. Although there is extensive 
literature describing hadronic transitions among bottomonium states nearly all of it is 
directed toward understanding the various tttt transitions. [2| While the rni transitions do 
provide important information about strong interaction dynamics in heavy quark systems, 
the investigation of other hadronic decay modes (i.e. involving rj, uo or multiple vr) should 
offer a different perspective. The analysis of heavy quarkonium hadronic transitions is one 
of few possible laboratories for the study of the physics of the soft gluon emission and 
hadronization process that governs such decays. In this letter, we describe the evidence for 
observation of new hadronic transitions, Xfei,2(2-P)— >co'T(lS). 

The data set analyzed for the purposes of this paper consists of (5.81 ± 0.12) x 10^T(3S) 
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decays observed with the CLEO III detector at the Cornell Electron Storage Ring. The 
CLEO III detector is described in Ref. Charged particle tracking is done by the 47-layer 
drift chamber and a four-layer silicon tracker which reside in a 1.5T solenoidal magnetic 
field. Photons are detected using an electromagnetic calorimeter (CC) consisting of 7784 
CsI(Tl) crystals distributed in a barrel geometry. The particle-identification capabilities of 
the CLEO-III detector are not used in the present analysis. 

The final state of 77r"'"7r~7r''£+£~ can be selected in a straightforward manner. Events with 
two high momentum charged tracks (p > 4 GeV) and two or three charged low momentum 
tracks (0.12 < p < 0.75 GeV) are considered. The low momentum tracks are required to 
come from the interaction region using criteria obtained by studying charged pion tracks in 
a sample of events from the kinematically similar decay T(2S)^7r"'"7r~T(lS). 

We require events to contain a T(1S) candidate by requiring that the two high momentum 
tracks in the event have an invariant mass consistent with the T(1S) mass (9.3 — 9.6) GeV. 
We make no cuts on the track quality variables for the lepton candidate tracks, nor do we 
require them to satisfy lepton-identification criteria. The invariant mass requirement alone 
provides a nearly background free sample, and imposition of lepton-identification criteria 
only leads to larger systematic errors and reduced signal efficiency without much improve- 
ment in signal quality. 

We require events to have three or four showers in the CC that are not matched to 
charged tracks. Two of these showers must form an invariant mass within 3a of the nominal 
tt" mass. In the further analysis, we utilize the momentum for this 7r° candidate obtained 
from a mass-constrained kinematic fit. In addition to the two showers that correspond to the 
tt'', events must contain an isolated photon candidate, between 50 and 250 MeV in energy, 
that does not form an invariant mass within 8 MeV ~ 1.5a of the vr*^ mass with any other 
shower. Further, the polar angle {6) of the third shower is restricted to | cos6'| < 0.804, the 
angular region in which CLEO's energy resolution is sufficiently good to resolve the photons 
for the signal channels of interest. Finally, we allow events to contain up to one additional 
shower in the range | cos^| < 0.804. 

In addition, we allow for the possibility of one spurious charged track candidate in ad- 
dition to the four "signal" tracks. Such spurious tracks may arise from failures in pattern 
recognition and from delta rays. Spurious showers may arise from synchrotron radiation 
from the beams or as a result of random noise in the CC. If a given event yields more 
than one candidate due to the presence of an additional shower or track, we choose the 
candidate for which the sum of energies of all final state particles is nearest M(T(3S)). 

Because there is no phase space for a pair of kaons for decays in which a T(1S) is present, 
we assume that the low momentum charged tracks are due to vr^. The invariant mass of 
the 7r+7r^7r° combination, plotted in Figure El shows a clear enhancement at the mass of u, 
M{uj) = 0.782 GeV. 

Finally, to complete the reconstruction of the full decay chain, 
T(3S)— s>7u;T(lS)-^77r+7r^7r°£+£~, we require that the be less than 10 for a kinematic fit 
of vr+vr^vr'^ constrained to the ut mass, and subsequently that the mass recoiling against the 
kinematically fitted u and 7 candidates lie within I20 MeV of M(T(1S)) = 9.460 GeV. 

The simplest explanation for the observed events is the decay sequence T(3S)^7Xbj(2P), 
with Xbj(2P)— *u;T(lS). Transitions through Xbi(2P) and Xb2(2P) only are allowed by energy 
conservation. xwC^P) lies below threshold for decay to co'T(lS). In principle a transition 
through the r]b{3S) state (see Figure H} is possible, but this state has never been observed, 
and, furthermore, the energy of the photon in the decay T{3S)^-yrih{3S) is expected to be 
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FIG. 2: vr+vr^vr invariant mass, for data events subject to final analysis cuts with the exception 
of the cut on the of the kinematic fit to u. The overlaid histogram shows signal Monte Carlo 
events (normalized to the same total number) and indicates the good reproduction by the Monte 
Carlo of the shape and location of the uj peak. 

below the range of observed energies in the data. We also neglect as highly improbable the 
possibility that the uj is emitted directly in a transition from T(3S) to a heretofore unknown 
state, which then decays radiatively to T(1S). 

Backgrounds from ordinary udsc quark production are extremely small because of the 
presence of the T(1S)— decay in the signal sample. The only significant source of 
background T(1S) expected is known cascades from T(3S). The final state of 77r+7r~7r°T(lS) 
may be reached through: 

T(3S) ^ 7X6(2P), X6(2P)^7T(2S), T(2S)^7r+7r-T(lS) 
T(3S) ^ 7rVT(2S), T(2S)^7r+7r-T(lS). 

In the first case, this final state can be produced by the addition of a spurious shower in the 
calorimeter. In the second case, it may be reached by loss of one of the 7 from one of the 
7r° due to acceptance or energy threshold. 

Backgrounds produced through cither of these two processes may be removed by excluding 
events in which the mass recoiling against the two charged pions in the T(3S) reference 
frame is consistent with the hypothesis that the tt^tt" system is recoiling against T(2S) in 
the process T(3S)^X + T{2S)^X + 7r+7r-T(lS) (between 9.78 and 9.81 GeV). 
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Two other decay sequences can yield the final state of 77r"'"7r 7r'^T(lS): the decay 
T(3S)^7r+7r-T(2S), with the T(2S) decaying either to 7r07rOT(lS), or to jXb{lP) followed 
by Xb{^P)^j'^ i^S) . In each of these cases, however, the vr^ are too soft to produce spurious 
u candidates for the signal decay chain. 

In order to evaluate residual background, we generated a GEANT Monte Carlo sample 
for the channel T(3S)^7Xb(2P), X6(2P)^7T(2S), T(2S)^7r+7r-T(lS), corresponding to 
21.5 ± 3.1 million T(3S) decays, or 4.53 ± 0.65 times our data set. The uncertainty on the 
equivalent number of T(3S) decays is due to the error on the branching ratios needed to 
convert our number of generated events to the equivalent number of T(3S) decays. This 
Monte Carlo sample produced one event that satisfied our selection criteria. We therefore 
expect 1/(4.53 ± 0.65) = 0.22 ± 0.03 events due to this source. We also generated a Monte 
Carfo sample of T(3S)-^7r07rOT(2S), T(2S)^7r+7r-T(lS), corresponding to 5381^?^ million 
T(3S) decays, or 113l;ig times our data set. From this sample, a total of nine events passed 
our selection. In our data set, we thus expect 9/(113 ± t^g) = 0.08 ± 0.01 events due to 
this source. To account for the residual background, we subtract the expected contribution 
of 0.30 events from the observed yield. We conservatively set a systematic error of ±0.30 
events due to this subtraction. 

To evaluate the signal detection efficiency, we generated 150,000 signal Monte 
Carlo events for each of Xfei(2-P) and X62(2P), proceeding through the sequence 
T(3S)— s>7Xw,2(2P)— >-7a;T(lS)— s>77r+7r~7r°£"''£~, and utilizing PDG average values for 
the X6i,2(2-P) and T(nS) masses, and with a uniform angular distributions for the 
T(3S)— >7Xbi,2(2-P) and T( IS)— decays. The analysis cuts described above have been 
applied to these samples, and we obtain e(xfe2(2-P)) = (6.23 ± 0.06)% and e(xhi(2P)) = 
(6.81 ± 0.07)%, including all selection criteria, acceptance and trigger efficiencies. The ef- 
ficiencies have been corrected for the known angular distributions of the El transitions 
T(3S)— s-7X6i,2(2P). We apply an additional relative systematic error of ^3% to the effi- 
ciency in order to account for the possibility that the T(1S) retains the initial polarization 
of the T(3S). 

In Figure El we present a scatter plot of the recoil mass against uj'-f vs. dilepton invariant 
mass for all events subject to all the cuts discussed above, except those on the variables 
plotted, in order to emphasize the clarity of the signal. The final spectrum is shown in 
Figure El The observed yield has possible contributions from both decay sequences involving 
X6i(2P) and X62(2P) intermediate states. 

To obtain branching fractions for the X62(2P) and X6i(2P) transitions, we perform a 
maximum likelihood fit of the E^ spectrum. The expected photon spectra for T(3S) transi- 
tions to Xfe2(2P) and X6i(2P) were obtained from the signal Monte Carlo samples, and the 
observed Ej spectrum was then fit to normalized Monte Carlo lineshapes with intensities 
(or yields) for XbiC^P) and Xfe2(2P) as the free parameters. We obtain yields of 32.6tg;i and 
20.1^5;! events, respectively. These yields have statistical significances of 10.2cr and 5.2a, 
respectively. The histogram resulting from the best fit is shown superimposed on the data 
in Figure El 

We subtract the expected background contribution of a total of 0.30 events, subtracting 
from each yield a fraction of this total background equal to the ratio of the fitted yield to 
the sum of the two yields. 

We thus obtain the following product branching ratios, using the detection efficiency and 
number of T(3S) discussed above: 

S(T(3S)-.7Xw(2P)) X B{xbi{2P)^iuT{lS)) 
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FIG. 3: Plot of recoil mass versus dilepton mass, for data events subject to the final set of cuts, 
with the exception of the cuts on the two variables plotted. 

xB{uj^n+7i-7T°) X i3(T(lS)^£+r) 

= 0.82+°;}^ ± 0.06 x 10-^ and (1) 
S(T(3S)^7X62(2P)) X i3(x62(2P)^c.T(lS)) 
xi3(cu^7r+7r-7r°) x i3(T(lS)^£+r ) 

= 0.55^[J;}^ ± 0.04 x 10~^ (2) 

in which the first error is statistical and the second is the systematic error. 

The statistical error of nearly 20% is the dominant source of error. Systematic error con- 
tributions to the above product of branching ratios are the following: 2% for the uncertainty 
in the number of T(3S), 1% per charged track (a total of 4%) for track finding, 5% for 7r° 
reconstruction, 2% for radiative 7 reconstruction, 1% for Monte Carlo statistics, to% for 
the assumption of uniform T(1S)— angular distribution, and 1% for background sub- 
traction. These contributions, added in quadrature, result in an overall relative systematic 

t +7.7 
error of Zj i- 

Using the present world average branching fractions for T(3S)— >7Xf)L2(2P), T{lS)^£~^i~ 
(taken to be twice that of T(1S)— and u^Tr'^n^n^ from Ref. jsf, we thus obtain 

S(x.i(2P)^c^T(lS)) = 1.63l°f and (3) 
B{x,2{2P)^ujT{1S)) = 1.10^°;i+°;}J%, (4) 

in which the errors quoted are statistical and systematic, respectively. The systematic 
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FIG. 4: Fitted photon energy spectrum for the final selection of events. The solid histogram shows 
contributions for both Xbii'^P) and Xb2(2-P), while the dotted and dashed histograms show the 
individual X6i(2-P) and Xf)2(2-P) contributions, repsectively. The dotted line indicates the region 
above which 7 energies are disallowed for lack of phase space. The small leakage of the histogram 
for the Xfei(2P) Monte Carlo sample into the kinematically disallowed region is due to the finite 
detector resolution. 



errors on include the experimental errors on the branching ratios for T(3S) ^7x51,2 (2P), 
u}^TT~^7r~7r^ and T(1S)— which contribute an additional systematic error of 5.9%. 

We may also calculate the ratio of Xb2{'2P) to Xbi(2P) branching ratios, for which sys- 
tematic errors discussed above cancel almost entirely. We obtain this through a maximum 
likelihood fit to the spectrum, in which the two free parameters are the sum of yields 
and the ratio of X62(2-P) to Xf,i(2P) yields. When this fit is performed, we obtain a sum of 
yields equal to 52.4tg;9 and a ratio of 0.62to;2o- order to convert the ratio of yields to the 
ratio of branching ratios, we multiply the yield ratio by a factor of e(x6i(2P))/e(xb2(2P)). 



B(x62(2P)^a;T(lS)))/B(x(.i(2P)^a;T(lS)) = 0.62tg 



e(x.2(2P)) 



= 0-62^S1Sx|f^^ = 0.67-30. 

(6.23 ±0.06)% ^ ^ 

The only systematic errors which do not cancel in this ratio are the small uncertainties in 
efficiency and T(3S)— >7X6i,2(2P) branching ratios. These are entirely negligible compared 
to the output error from the maximum likelihood fit. 
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In Ref. p], Voloshin predicts on the basis of S-wave phase space factors, that the ratio 
of partial widths r(x;,2(2P)— i^a;T(lS))/r(x;,i(2P)— i>ci;T(lS)) should be approximately 1.4. 
The ratio of full widths r(xb2(2P))/r(xfci(2P)) lies in the range of 1/(1.5-1.25), using world 
average measurements of i3(xfei^2(2-P)— >7T(1S, 2S)) and theoretical predictions for the rates 
r(x;,i,2(2P)^7T(lS, 2S)). Thus the branching ratios S(x6i,2(2-P)— >c<jT(1S)) are expected 
to be approximately equal. Our measurement is rough agreement with this expectation. 

Furthermore, Gottfried i3] predicted that, due to the non-Abelian nature of strong in- 
teractions, AL = 1, AS' = hadronic transitions ought to have sizeable widths. Our 
observation provides experimental evidence for his statement. 

We have made the observation of the first hadronic decay modes involving non-vector bb 
states, thanks to the large statistical advantage of the CLEO-III bb data set over all previous 
data sets. The discovery came in the phase-space limited channels Xfei,2(2-P)— >cijT(1S), and 
the ratio of the measured branching ratios for the two transitions are in rough agreement 
with theoretical expectations based on S-wave phase space factors. The observation of these 
decays supports the prediction that such transitions should have significant rates. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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